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Abstract 

Coherent motion of cold atoms in a standing-wave field is interpreted as 
a propagation in two optical potentials. It is shown that the wave-packet 
dynamics can be either regular or chaotic with transitions between these 
potentials after passing field nodes. Manifestations of de Broglie-wave chaos 
are found in the behavior of the momentum and position probabilities and 
the Wigner function. The probability of those transitions depends on the 
ratio of the squared detuning to the Doppler shift and is large in that range 
of the parameters where the classical motion is shown to be chaotic. 
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1. Introduction 

Cold atoms in a standing-wave laser field are ideal objects 

for studying fundamental principles of quantum physics, quantum-classical 
correspondence, and quantum chaos. The proposal [f| to study atomic dy- 
namics in a far-detuned modulated standing wave made atomic optics a test- 
ing ground for quantum chaos. A number of impressive experiments have 
been carried out in accordance with this proposal 0, S, @]. Furthermore, 
cold atoms have been recently used to study different phenomena in statis- 
tical physics including ratchet effect of a directed transport of atoms in the 
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absence of a net bias I 111 . |12| and Brillouin-like propagation modes in 
optical lattices [13 



New possibilities are opened if one works near the atom-field resonance 
where the interaction between the internal and external atomic degrees of 
freedom is intense. A number of nonlinear dynamical effects have been found 



141 . Il5l . [16| with point-like atoms in a near-resonant rigid optical lattice: 
chaotic Rabi oscillations, chaotic walking, dynamical fractals, Levy flights, 
and anomalous diffusion. 

Dynamical chaos in classical mechanics is a special kind of random-like 
motion without any noise and/or random parameters. It is characterized by 
exponential sensitivity of trajectories in the phase space to small variations 
in initial conditions and/or control parameters. Such sensitivity does not 
exist in isolated quantum systems because their evolution is unitary, and 
there is no well-defined notion of a quantum trajectory. Thus, there is a 
fundamental problem of emergence of classical dynamical chaos from more 
profound quantum mechanics which is known as quantum chaos problem 
and the related problem of quantum-classical correspondence. In a more 
general context it is a problem of wave chaos. It is clear now that quantum 
chaos, microwave, optical, and acoustic chaos have much in common (see 



17( for a review). The common practice is to construct an analogue for a 
given wave object in a semiclassical (ray) approximation and study its chaotic 
properties (if any) by well-known methods of dynamical system theory. Then, 
it is necessary to solve the corresponding linear wave equation in order to 
find manifestations of classical chaos in the wave-field evolution in the same 
range of the control parameters. If one succeeds in that a quantum-classical 
or wave-ray correspondence is announced to be established. 

In this paper we perform this program with cold two-level atoms in a 
one- dimensional standing-wave laser field and show that coherent dynamics 
of the atomic matter waves is really complicated in that range of the con- 
trol parameters where the corresponding classical point-like atomic motion 
can be strictly characterized as a chaotic one. The effect is explained by a 
proliferation of atomic wave packets at the nodes of the standing wave. 



2. Order and chaos in dynamics of atomic wave packets in a laser 
standing wave 

The Hamiltonian of a two-level atom, moving in a one-dimensional clas- 
sical standing-wave laser field, can be written in the frame rotating with the 



2 



laser frequency ojf as follows: 



P 2 h 

H = - h ~{oJ a - Uf)<J z - Ml (<x_ + a + ) cos kfX, 

—,11 Lfi 



(1) 



where a± tZ are the Pauli operators for the internal atomic degrees of freedom, 
X and P are the atomic position and momentum operators, u a and Q are 
the atomic transition and Rabi frequencies, respectively. We will work in the 
momentum representation and expand the state vector as follows: 



*(«)>= [a(P,t)\2) + b(P,t)\l)]\P)dP, 



(2) 



where a(P, t) and b(P, t) are the probability amplitudes to find atom at time 
t with the momentum P in the excited and ground states, respectively. The 
Schrodinger equation for these amplitudes is 



id(p) = X -{uJrV 2 ~ A)a(p) - -[b(p + 1) + 6(p - 1)], 
ib(p) = ^{u rP 2 + A)6(p) - -[a(p + 1) + a(p - 1)], 



(3) 



where dot denotes differentiation with respect to dimensionless time r = fit, 
and the atomic momentum p is measured in units of the photon momentum 
Hkf. The normalized recoil frequency, u r = hk 2 j-/m a Q, and the atom-field 
detuning, A = (cuf — uj a )/VL, are the control parameters. 



We will treat the wave-packet motion in the dressed-state basis [18 

| + ) A = sin6|2) + cos6|l), | - ) A = cos 0|2> - sin 6|1). (4) 

The dressed states are eigenstates of an atom at rest in a laser field with the 
eigenvalues of the quasienergy and the mixing angle 0: 
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(5) 



The depth of the nonresonant optical potential is 

U A = 



|A| 



(6) 
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At A = 0, the depth changes abruptly from 1 to 2. The resonant 
(A = 0) and non-resonant (A ^ 0) potentials of an atom in a standing- 
wave field are drawn in Fig. [TJ The ground atomic state can be written as a 
superposition of the dressed states 

ii\ i \ 1 i \ tan 6 . . 

V 1 + tarn v 1 + tarn 9 

In general, an atom in the ground state, placed initially at xq = 0, will move 
along two trajectories simultaneously because it is situated simultaneously 
at the top of E^ and the bottom of E^ (see FigJT]). 

In the dressed-state basis, the probability amplitudes to find the atom at 
the point x in the potentials E^ and E^ ^ are, respectively, 

C+(x) = a(x) sin 6 + b(x) cos 6, C_(x) = a(x) cos 6 — b(x) sin 6, (8) 

where the amplitudes in the bare-state basis, a(x) and b(x), are computed in 
the position representation with the help of the Fourier transform 

/OO /'OO 
dp'e ip ' x a(p'), b(x) = const / dp' e ip ' x b{p'). (9) 
-OO J — OO 

To clarify the character of the bipotential motion we write down the 
Hamiltonian of the internal degrees of freedom in the basis | ± ) = ( 1 1) ± 
|2»/V2: 

A . . 

Hi n t = cr z cos x + —a x . (10) 

Let us linearize the potential in the vicinity of a node of the standing wave 
and estimate a small distance the atom makes when crossing the node as 
follows: 5x ~ w r bnodc| r , where |p no de| is the mean atomic momentum near 
the node. The quantity w r |p no d e | is a normalized Doppler shift for an atom 
moving with the momentum (pnodel, i-e., (jJd = ^ri^nodel = kf\v no & e \/Vt. We 



arrive now at the famous Landau-Zener problem [19] to find a probability of 
transition between the states (j4j) when the energy difference varies linearly 
in time. In other words it is a a probability of transition from one potential 
to another or from one trajectory of motion to another. The asymptotic 
solution is 

ttA 2 ' 



P LZ = exp( ). (11) 

Up 



There are three possibilities. 
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1. A 2 ^> u>d- The transition probability is exponentially small and an 
atom moves adiabatically along, in general, two trajectories without any 
transitions between them. 

2. A 2 <C cj£). The probability of a Landau-Zener transition is close to 
unity, and an atom changes the potential each time upon crossing any node, 
i.e., it moves mainly in the resonant potentials E^. 

3. A 2 ~ uj£>. The probabilities to change the potential or to remain in 
the same one upon crossing a node are of the same order. In this regime one 
may expect strong complexification of the wave function. 

At large and small detunings, the translational motion splits into two 
independent motions in the potentials and the wave-packet motion is 

regular in the first two cases. In the third case, the motion is complex because 
of a proliferation of wave packets at the nodes of the standing wave. We call 
such a motion as " a chaotic" one by the reasons which will be clear in section 
3. The switch between the regular and chaotic regimes of atomic motion can 
be easily performed by changing the detuning. 

Our normalization enables to change dimensionless values of the recoil 
frequency u r and the detuning A varying the Rabi frequency Q. Working 
with a cesium atom at the transition 6S1/2 - 6P3/2 (m> a = 133 a. u., A a = 852.1 



nm and z/ rec — 2 KHz 20j), we have u r = 10 3 at Q = 1 MHz. Let the atom 



is initially prepared in the ground state as a Gaussian wave packet in the 
momentum space with p = 55. The probability to find the atom with the 
momentum p at time r is P(jp, r) = \a(p, r)| 2 + \b(p, t)| 2 . 

To illustrate the difference between the regular and chaotic regimes of 
the wave-packet motion we take two values of the detuning, A = 1 and A = 
0.2. At A = 1, the motion is expected to be adiabatic and regular because 
A 2 u D = 0.055, and the Landau-Zener probability, P LZ , is exponentially 
small. At A = 0.2, one expects a much more complicated wave-packet motion 
with nonadiabatic transitions between the two potentials at the field nodes 
because A 2 = 0.04 ~ ojd- Figure |2] shows the dependence of the mean atomic 
momentum < p > over a large time scale in those two cases. In both the 
cases < p > oscillates in a rather irregular way. The difference is that for an 
adiabatically moving wave packet, which we refer as a regular motion, the 
mean atomic momentum oscillates in a narrow range around po = 55 (the 
upper curve in the figure). Whereas the range of its oscillations for a wave 
packet moving with nonadiabatic transitions at the field nodes, which we 
refer as a chaotic motion, is much more broad (the lower curve in the figure). 
It is a simple illustration of the two different regimes of the wave-packet 
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propagation in terms of the classical variable. 

In Fig. [3k we plot the dependence of the momentum probability-density 
on time at A = 1. The initial wave packet splits from the beginning to a few 
components because the initial ground state is a superposition of the dressed 
states (j3J). The initial kinetic energy is enough to perform a ballistic motion. 
The momentum changes in a comparatively small range, from 40 to 70 of 
the photon-momentum units. The packet does not split at the nodes of the 
standing wave but it, on the contrary, recollects in the momentum space at 
the nodes and spreads in between. However, this recollection smears out in 
course of time. 

At A = 0.2, the atomic ground state is the following superposition of 
the dressed states: |1) ~ 0.74| + ) A + 0.66 1 — ) A . The | + ) A -component of 
the initial wave packet, i.e., that one, starting from the top of the potential 
overcomes the barriers of that potential and moves in the positive di- 
rection of the axis x proliferating at the nodes. As to the | — )-component 
with decreased values of p, it will be trapped in the potential well perform- 
ing oscillations in the momentum and position spaces. The period of those 
oscillations is about T ~ 280 which is equal approximately to the period 
of revivals of the Rabi oscillations for the population inversion. Figure [3b 
illustrates the effect of simultaneous trapping and ballistic motion of the 
atomic wave packet in the chaotic regime resulting in a broad momentum 
distribution, from p = —60 to p ~ 80. 

To illustrate the nonadiabatic transitions from one potential to another 
and their absence at the nodes more explicitly, we go to the position space 
and compute the probabilities \C±(x,r)\ 2 (jSJ) to be at the point x at time 
t in the potentials and £7 A \ respectively. In Fig. [4] we plot the evo- 
lution of the probability density |C_(x,r)| 2 in the frame moving with the 
initial atomic velocity oj r Po — 0.055. The slope straight lines in the figure 
mark positions of the nodes in the moving frame. In the case of the regular 
motion at A = 1 (Fig. H^), no transitions happen when the atom crosses 
the nodes. In the chaotic regime at A = 0.2, one observes visible changes 
in the probability-density |C_(x)| 2 exactly at the node lines (see Fig. 0b). 
It means transitions from one trajectory to another at the field nodes that 
should occur in a specific range of the control parameters if A 2 ~ up. This 
results in a proliferation of the components of the wave packet at the nodes 
and, therefore, a complexification of the wave function (see Fig. |3b). 

The Wigner function can be used to visualize complexity of the wave func- 
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tion in the chaotic regime of the atomic motion. We compute the evolution 
of the Wigner function of the ground state in the momentum space 

/oo 
dp'e' ip ' x a(p - p'/2)a*(p + p'/2). (12) 
-oo 

This quantity gives a quasi-probability distribution corresponding to a gen- 
eral quantum state ([2]). Figure [5] shows a contour plot of the Wigner function 
f[T2"j) at two moments of time, r = 50 and r = 200, when the atom moves in a 
regular way (A = 1). Figure [6] is a contour plot of this function at the same 
times, but with an atom making nonadiabatic transitions at the field nodes 
(A = 0.2). In the chaotic regime (Fig. [6Jd) we see a dust-like distribution 
of nonzero values of the Wigner function at r = 200 which occupy much 
more larger area in the phase space than the function for the regular motion 
(Fig. [5b). 



3. Quantum-classical correspondence 



In this section we compare the results of the quantum treatment with 
those obtained for the same problem in the semiclassical approximation when 



the translational motion has been treated as a classical one [14], [15j. We must 
compare quantum results for a single atomic wave packet not with a single 
point-like atom but with an ensemble of point-like atoms. Dynamical chaos 



has been found and analyzed in detail in Refs. [14j, [15j in the semiclassical 
approximation. Both the internal and external degrees of freedom of a two- 
level atom in a standing wave field have been shown to be chaotic in a specific 
range of values of the detuning A, the recoil frequency u r , and the initial 
momentum p . 

Coherent semiclassical evolution of a point-like two-level atom is governed 
by the Hamilton-Schrodinger equations with the same normalization as in the 
quantum case 



15 



x = (jJ r Pi P=— usinx, u = Av, 
) = — Am + 2z cos x, z — —2v cos x, 



(13) 



where 

u = 2Re(a 6*), v = -2lm(a b* ), z = |a | 2 - |6 | 2 (14) 

are the atomic-dipole components (u and v) and population-inversion (z), 
and ao and 6q are the complex-valued probability amplitudes to find the 
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atom in the excited, |2), and ground, |1), states, respectively. The system 
( II 3p has two integrals of motion, the total energy 



and the length of the Bloch vector, u 2 +v 2 +z 2 = 1. Equations ([TBI) constitute 
a nonlinear Hamiltonian autonomous system with two and half degrees of 
freedom and two integrals of motion. 

In classical mechanics there is a qualitative criterion of dynamical chaos, 
the maximal Lyapunov exponent A, which measures the mean rate of the 
exponential divergence of initially closed trajectories in the phase space. In 
Fig. [7]we plot this quantity, computed with semiclassical equations of motion 
(fl3"j) vs the detuning A at w r = 1CT 3 and po = 55. At zero detuning, the 
set of semiclassical equations acquires an additional integral of motion and 
becomes integrable. The center-of-mass motion is regular at small values of 
the detuning, A C 1, and at large ones, A > 0.8. Positive values of A at 
< A < 0.8 characterize unstable motion. Local instability of the center- 
of-mass motion produce chaotic motion of an atom in a rigid standing wave 
without any modulation of its parameters in difference from the situation 
with atoms in a periodically kicked optical lattice 0, There is a range 
of initial conditions and the control parameters where the center-of-mass 
motion in an absolutely deterministic standing wave resembles a random 
walking. It means that a point-like atom alternates between flying through 
the lattice, and being trapped in its wells changing the direction of motion 



In Fig. [8] we illustrate semiclassical chaos with the Poincare mapping 
for a number of atomic ballistic trajectories in the western (u < 0) and 
eastern (u > 0) hemispheres of the Bloch sphere (u, v, z) on the plane v — z. 
One can see a typical structure with regions of regular motion in the form 
of islands and chains of islands filled by regular trajectories. The islands 
are imbedded into a stochastic sea, and they are produced by nonlinear 
resonances of different orders. Increasing the resolution of the mapping, one 
can see that large islands are surrounded by islands of a smaller size each of 
which, in turn, is surrounded by a chain of even more smaller islands, and so 
on. 

In quantum mechanics there is no well-defined notion of a trajectory in 
the phase space, the very phase space is not continuous due to the Heisen- 
berg uncertainty relation, and, hence, the Lyapunov exponents can not be 



W = — p —ucosx z, 



(15) 



in a random-like way [15|, 
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computed (however, see Ref. [21( where a notion of a Lyapunov exponent for 
quantum dynamics has been discussed). The main result of this paper is the 
establishment of the fact, that chaotically-like complexification of the wave 
function, caused by nonadiabatic transitions at the field nodes, occurs exactly 
at the same range of the control parameters where the semiclassical dynamics 



has been shown to be chaotic in Refs. [14J, |15(. It should be stressed that 
quantum motion of a wave packet with nonadiabatic transitions between the 
two optical potentials is compared with the center-of-mass motion of an en- 
semble of atoms each of which moves in a single optical potentials. So, when 
we say about a quantum-classical correspondence we mean a correspondence 
between the wave function of a single quantum atom and the trajectories of 
the ensemble of classical atoms with different values of the initial momentum 
Po and at the other equal conditions. 

4. Conclusion 

We have studied coherent dynamics of cold atomic wave packets in a 
one-dimensional standing-wave laser field. The problem has been considered 
in the momentum representation and in the dressed-state basis where the 
motion of a two-level atom was interpreted as a propagation in two optical 
potentials. The character of that motion has been shown to depend strongly 
on the ratio of the squared detuning, A 2 , to the normalized Doppler shift 
u>d- In the regular regime, when A is comparatively large or small, wave 
packets move in a simple way. The chaotic regime occurs if A 2 ~ wo when 
the probability for an atom to make nonadiabatic transitions while crossing 
the nodes of the standing wave is large. Atom in this regime of motion 
simultaneously moves ballistically and is trapped in a well of the optical 
potential. This type of motion and proliferation of wave packets at the nodes 
result in a complexification of the wave function both in the momentum and 
position spaces manifesting itself in the irregular behavior of the Wigner 
function. 

Comparing the results of the quantum treatment with those obtained 
in the semiclassical approximation, when the translational motion has been 



treated as a classical one [14], [15] , we have found that the wave-packet dy- 
namics is complicated exactly in that range of the atom-field detuning and 
recoil frequency where the classical center-of-mass motion has been shown to 
be chaotic in the sense of exponential sensitivity to small variations in initial 
conditions or parameters. 
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As to possible practical applications of the results obtained we mention 
atomic lithography to produce small-scale complex prints of cold atoms (see, 
for example, beautiful experiments on coherent matter-wave manipulation 



22l . l23l . l24j), new ways to manipulate atomic motion in optical lattices by 



varying the atom-filed detuning and atomic ratchets with cold atoms. 
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Figure 1: The resonant E^ (A = 0, dotted curves) and non-resonant E^ (A 7^ 0, 
solid curves) potentials of a two-level atom in a standing-wave laser field. An atomic wave 
packet, centered at xo = 0, and its initial evolution in the upper and lower potentials are 
shown schematically. 




Figure 2: Mean atomic momentum (p) vs time in the regular (A = 1, the upper bold 
curve) and chaotic (A = 0.2) regimes of motion. 
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Figure 3: Momentum probability-density distribution P(p,r) vs time in (a) the regular 
and (b) chaotic regimes of motion. The color codes the corresponding values of P(p). 




Figure 4: The probability |C_(a;)| 2 to find the atom in the potential E A in the moving 
frame of reference in (a) the regular and (b) chaotic regimes of motion. The slopes mark 
positions of the nodes in the moving frame and the color codes the values of \C-(x)\ 2 . 
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Figure 5: Contour plots of the Wigner function at (a) r = 50 and (b) r = 200 for the 
regular regime of the atomic motion. Color online: red and blue areas show positive and 
negative values of the Wigner function, respectively. 
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Figure 6: The same as in Fig. [5] for the chaotic regime of motion. 
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Figure 7: Maximal Lyapunov exponent A, computed with semiclassical equations of motion 
(fT3|) . vs the detuning A at u) r = 10~ 3 and po = 55. 
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Figure 8: Poincare sections of the Bloch sphere illustrating the effect of semiclassical chaos 
with point-like atoms at u> r = 10~ 5 , A = —0.05 and the total atomic energy W = 36.45. 
(a) u < (western Bloch hemisphere), (b) u > (eastern Bloch hemisphere). 
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